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I We investigate the phase diagram of a quarter filled Hubbard ladder with nearest-neighbor 

Coulomb repulsion using bosonization and renormalization group approach. Focusing on the strong- 
' repulsion regime, we discuss the effect of an interchain exchange interaction J± and interchain re- 

' pulsion V± on the possible ground states of the system and charge order configurations. Since the 

, spin excitations always possess a gap, we find competing bond-order wave and charge density wave 

H ' phases as possible ground states of the ladder model. We discuss the elementary excitations in these 

various phases and point an analogy between the excitations on some of these phases and those of 
5^ , a Kondo-Heisenberg insulator. We also study the order of the quantum phase transitions between 

■ the different ground states of the system. We obtain second order transitions in the Ising or 5(7(2)2 

^ ; ' universality class or first order transitions. We map the complete phase diagram in the J± — V± 

, plane by integrating perturbative renormalization-group equations. Finally, we discuss the effect of 

. doping away from half-filling and the effect of an applied magnetic field. 

^ : 

P3 . I. INTRODUCTION 

o : 

1 ^ I ' Charge ordering is a general phenomenon in condensed matter physics that has been recently observed in a variety of 

I compounds including rare earth manganitesi*^, quasi-one-dimensional (TMTTF)2X conductors^, cuprate or nickelate 
I materials^. A very particular form of charge order in these latter materials is the formation of charge stripes, that 
^ I are domain walls between hole rich and hole poor regions. From the theoretical point of view, charge ordering is a 
. crystallization of the electron liquid which occurs when long-range Coulomb repulsion^ dominates over kinetic energy. 
' In one dimension, charge ordering and more generally metal-insulator transitions can be studied in great details due 
, to the existence of powerful analytical2i2iS and numericaliS methods. In the case of a one-dimensional system, the 
' simplest model^^ of interacting electrons that allows for charge ordering is the Hubbard model extended with an 
additional nearest-neighbor Coulomb interaction V. Studies of the metal insulator transition in one dimension show 
that charge ordered or Mott insulating states can form at commensurate fillings for sufficiently strong repulsioniSiH. 
In relation with organic compounds, the quarter- filled extended Hubbard model on a single chain has been studied in 
detailsii. As a first step towards understanding charge ordering and stripe formation in strongly correlated electron 
Ch ' models, numerical studies of coupled chain systems at commensurate fiHingsiSiiSiiLiSiiS (ladders) have also been 
, performed. Half filled ladder systems are Mott insulators that display a spin gap analogous to the Haldane gap in 
(— I ' spin-1 chains^^ when made of an even number of coupled chains. Ladder models are not only of theoretical interest. 
Q Quite a few half-filled ladder systems have been synthesized both organic and inorganic and the spin gap behavior 
O ' has been characterized in great details^iiS2iSi24i2SiSSi. Away from half-filling and more generally commensurate filling, 
^ ] ladders are expected to be conducting. The persistence of the spin gap away from half-filling has been proposed 
to give rise to a paring mechanism and to superconductivity^S. This suggestion has given rise to intense theoretical 
studies both analyticali2i224^2i^ and numerica]p224Si24i2Si2& that confirmed the presence of superconducting correlations 
^ ■ at incommensurate fillings. Experimentally, the two-leg ladder compound Sri4Cu2404i— can be doped away from 
- - ' half filling and under pressure displays a superconducting transition. Recently, the compound a' — NaV2 05, initially 
identified as an inorganic spin Peierls compound^^, was shown to be a quarter- filled ladderi^SiiSiiiiiS, that undergoes 
a transition at Tc — 3AK corresponding to the formation of a charge ordered state with a zig-zag charge pattern and 
the opening of a spin gap^^i^^*^^*^. This has prompted numerical studies of quarter-filled extended Hubbard two-leg 
ladder systemaiL^ that exhibit the formation of a zig-zag charge ordered state at large repulsion. Interestingly, 
the charge ordered state of the two-leg ladder at quarter filling presents a spin gap in contrast to the single chain. 
Analytic studies of insulating states on a two-leg ladder have been mostly confined to the half-filled case^2i^2*^ii^ 
or to mean field^SiS approximations and strong coupling expansionsMi^^iS in the quarter-filled case. In the present 
paper we analyze the phase diagram of the quarter-filled two-leg Hubbard ladder with interchain Coulomb repulsion 
V± and exchange coupling J±, by using bosonization and renormalization group (RG) methods. Focusing on the 
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regime of strong repulsion, we discuss charge orders and the elementary excitations of the insulating phase. Due to 
the presence of a spin gap, the competing ground states are bond-order wave states (BOW) and charge density wave 
states (CDW). We analyze the phase transitions between these states using a perturbative renormalization group 
approach and a mapping of spin excitations to a Majorana fermion theory. The plan of the paper is the following. In 
Sec. m we introduce the model and discuss the physics of two simple limits, the large on-site repulsion and the large 
interchain exchange limit. In Sec. IIIII we give the details of the bosonization treatment and of the derivation of the 
renormalization group (RG) equations. In Sec. IIVI we introduce the order parameters of the insulating phases and 
the describe the charge ordered ground states obtained by the values of the phase fields that minimize the energy. 
In Sec. we describe the elementary charge and spin excitations of each ground state and discuss the analogy with 
Kondo-Heisenberg-Hubbard chain in a particular case. In Sec. I VII we discuss the nature of the transitions among the 
various phases in the antisymmetric charge sector and the spin sector. This is partially accomplished by a mapping to 
a theory of Majorana fermions and by a mapping to an effective quantum Ising model in the limit of strong interchain 
repulsion. Finally, we discuss the results of the phase diagram obtained by the numerical integration of the RG 
equations. We also briefly discuss the commensurate-incommensurate transitions produced by the application of a 
magnetic field or by doping away from quarter-filling. 

II. HAMILTONIAN AND SOME SIMPLE LIMITS 

We consider the quarter-filled extended Hubbard model on a two leg ladder. The Hamiltonian reads: 
H = -^5I("^J+i,p,'^'=''P'<^+4p,'TC«-Hi,P,'T)-i-LXl(42,^c,,i,<,-^H.c.) 

i,p i,p i 

where rii^p^a = clp o-Ci.p.cr, Tii^p = ?^i,p,| -I- ^i^p^i, i is the site index, p the chain index, t the intrachain hopping, t± the 
interchain hopping, U the on-site repulsion, V\\ the nearest-neighbor repulsion, and V± the interchain repulsion. 

Since we will be focusing on the case of very strongly repulsive interactions, the single-particle interchain hopping 
term t± will be irrelevant in this regime. However, it is well known that even when t± is irrelevant, it generates 
an interchain exchange term in the Hamiltonian22i^. Thus, the Hamiltonian should also include an interchain 
exchange term, 

Hexch. = J± Si,! • Si^2, (2) 
i 

where Si^p = c\ p ^da^pCi^p^jj. In the limit of t±^ ^ U, J± is estimated to second order perturbation theory in t±/U 
as Jj_ ~ t\/U ■ In the following, we will consider the Hamiltonian Q completed by the term Q and neglect 
altogether the interchain hopping t±. In fact, another model of interest is the two leg t-J ladder at quarter filling. 
The Hamiltonian of this model reads: 

H = -t^^ 'P{4+l.,p.aCi,p.a + c\p,^Ci+i^p^„)V - XI ^(42,aC,,i,<, + YL.C.)V 

i,p,cr i.a 

'^^,l?^^,2, (3) 

i,p i i 

where V is an operator that projects onto singly occupied statesi2i^. We will thus treat J± as a parameter independent 
of ?7 in (QJ in order to be able to discuss also results relevant to the two-leg t — J ladder model. In the rest of the 
present section, we will illustrate two simple limiting cases of the problem defined by the Hamiltonian (|2Jl which 
display insulating charge ordered ground states. 

A. large on-site repulsion 

Let us begin by considering the limit i7 ^ oo in the Hamiltonian ©-lEl- In this limit, it is not possible to put two 
fermions on the same site even when they have opposite spins and from the point of view of charge excitations, the 
system behaves as if it was made of spinless fermions with a Fermi wavevector kp twice the one of the spinful fermions 
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at [/ = 0^. The original spinful fermions system being quarter-filled, the spinless fermions system is half-filled and 
its effective Hamiltonian reads: 



(4) 



where now rii^p = alpUi^p. The bosonization of the Hamiltonian Q is straightforward. The relevant formulas can be 
found for instance in^. We obtain the following bosonized Hamiltonian: 



r=l,2 



dx 
2^ 



2g 



(27ra)2 



dx cos ^ 



dx cos 2(j)i cos 2(/)2 



Via 



dxaj;(/>i5^02,(5) 



where cos 401^2 represent the intrachain Umklapp interactions. Let us first consider the case in which these processes 
are irrelevant {V± — 0). This case corresponds to |V||| < 2t. The Hamiltonian Q has then a Luttinger liquid ground 
state for V±_ = 0, and the renormalized Luttinger exponent K* can be obtained from studies of the t — V model^ as: 



K* 



2-4: arccos ( 



(6) 



In the Luttinger liquid regime, we can set g — and replace K by K* in the Hamiltonian jSJ provided V± is sufficiently 
small. The Hamiltonian |(SJ, is then decoupled by introducing the fields (j)± — {<j)i ± <j)2)/V2 leading to the following 
Hamiltonian: 



H = H++H^, 



Hr = 



where r = ± and: 



dx 
2^ 



— u 1 + r 



2V^a 
(27ra)2 



KV^a 



dx cos 



(7) 



(8) 



For Kr < 1, the term cos vSi/)^ is relevant and opens a gap in the charge modes. If VI <C u/a, this implies that the 
gap opens as soon as K* < 1, whereas the intrachain Umklapp processes become relevant only for K* = 1/2. We 
thus see that for a wide range of < V|| < 2t, although intrachain repulsion V|| alone is too weak to open a charge 
gap by itself, the existence of a nonzero interchain repulsion is sufficient to induce an insulating ground state. This 
is consistent with the numerical observation of a charge gap state forming for V|| > and t± < t a,t large U in Rcf.^^. 

In that insulating state, the charge gap should vary as Ap+ ~ u/a{V±a/uy/''^~^+\ 

We note from ((HJ that K* = 1 corresponds to V\\ ~ 0, i.e. the Hamiltonian Q is identical to the bosonized 
Hamiltonian of the half- filled ID Hubbard chain, (j>^ playing the role of (j)p, and 0_ playing the role of (pa- As a 
result, although the total charge mode is gapful, the antisymmetric charge mode remains gapless. This insulating 
state presents the same SU{2) x SU{2) ^ 80(4) symmetry as the Hubbard model. In particular, the gapped charge 
excitations behave like those of the half-filled Hubbard model and the charge gap is Ap_|_ ~ u/aexp{—Cu/{V±a)). 

To determine the long range order that is realized in the ground state, we need to fix the values of {(t>i^2) that 
minimize the classical ground state energy, i. e. require that cos2((/)i) cos2((/)2) is negative. Since the density of the 
spinless fermions reads: 



Pp{x) — dx4>p{x) -I cos(20p — 2kpx) 



1 

na 



(9) 



where p = 1,2 is the chain index, this leads to two out of phase charge density waves of wavevector 2k'p = - on 
chains 1 and 2, i.e. a zig-zag charge ordering. This ordering is represented on figure ^ We note that in Ref. l48l it 
was found that the charge ordering was formed for U = oo only when V|| > 2t although a charge gap was obtained for 
1^1 > 0. The reason for this discrepancy in Ref. li^ could be the presence oit/tj_ ~ 1. Turning to transport properties, 
from lO the long wavelength density on chain p is Pp{x) = —-^dxipp- By considering the topological charge of the 
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sine-Gordon models one easily obtains that a charge solitons carries the electrical charge ±e. The a.c. electrical 
conductivity of this system can be obtained from the form factor expansionffiiSi. 

When Vji > 2t, the interchain Umklapp processes become relevant and induce the formation of charge density waves 
on each chain^^. The effect of V± is to lock the respective phases of these charge density waves into a zig-zag pattern. 
Zig-zag charge ordering is thus ubiquitous in the limit U = oo. 



B. large interchain exchange 

In the limit in which J± —> oo in ^ or in (PJ, it is again possible to derive a simplified Hamiltonian in a 
low-energy subspace. Namely, one can restrict to a subspace in which fermions form spin singlet pairs on the same 
rung. These pairs can then be treated as hardcore bosonaM moving on a single chain and carrying a charge 2e. For 
a quarter filled fermion system, the effective boson system is at half-filling. The effective bosonic Hamiltonian reads: 

H = -iJ2ibh+i+tl- c.)-Ky^n,n,+i, (10) 

i i 

where 6j = '^„o'Ci^i^aCi^2,-a-, annihilates a singlet pair on site i. In the limit J± ^ t, the coefficients t,V can be 
obtain from second order perturbation theory as i = t"^ / J±, V — 21"^ / J± + V\\. The Hamiltonian (fTUIl is solved exactly 
by the Bethe AnsataSSi, and its spectrum has a gap for Vj| > 0. A bosonization descriptionS2iSiS4 of the low energy 
excitations of this system gives: 



H 



dx 
2^ 



2A 

{2-Kay 



dxcosA4>b, (11) 



where the mode describes the charge excitations of the hardcore boson system. The electrical charge density is 
given by: 

2e e^'^^ 

p(x) = 9:r06 + 2e cos2(/>h (12) 

TT 7ra 

In the gapped regime, (0^) — 0, and the ground state of this system is formed of singlet pairs in a charge ordered 
state as represented on Fig. [7| We therefore see that in the regime of strong interaction, both Jl and will induce 
a charge ordering. However, Jl favors a phase with stripes formed along the rungs of the system, whereas Vj_ favors 
a zig-zag charge ordering. As a result, the competition between J_l and Vj_ induces a frustration in the system that 
can lead to a variety of charge ordering patterns. 



III. BOSONIZATION DESCRIPTION 



In the present section, we derive a bosonized representation of the Hamiltonian (E3)~(|2Jl in the limit Jl, Vj^ ^ t- For 
J_L = Vj^ = 0, the chains arc decoupled, and the Hamiltonian describing their low energy, long wavelength excitations 
reads: 



Hp,p ^ I dx 



dx 



H = ^ Hp^p + iJo-,j 

p=l,2 

UpKp{-KWp^pf + ^{dxcf)p.pf 



(13) 

(14) 
(15) 



where the fields satisfy to canonical commutation relations [(l>i,^p{x), Tl^'^p' (x')] = i5vy5p,pi5{x — x'), {v = p,a). Since 
we do not assume U,Vh <C t, we use the renormalized values of Up^ Kp^u^, in the bosonized Hamiltonian of the 
decoupled chains (|14|l - l|15() . Spin rotational invariance imposes the renormalized value of K^^ = 1. The renormalized 
values of the remaining quantities can be obtained from the Bethe Ansatz^^'SSiSLS for the Hubbard model or the t-J 
model at the supersymmetric point, or from numerical calculations in the case of a non-integrable model22iZ2iii. 
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A. derivation of interchain coupling 

In Sec. Ill Al we have seen that the Akp harmonics in the fcrmion density play a cr ucial ro le for U/t ^ 1. The 
expression of the Akp harmonics in terms of the boson fields can be obtained from Refs. 



Pr{x) 



V2 



+ 



(16) 



7ra 7ra 

where r = 1, 2 and kp — j^. For U —> oo, we recover the bosonized expression ^ of the density of spinless fermions 
with k'p = 2kp, (/) — \/20p and K = 2Kp. Using the expression IjlGf) of the fermion density, we obtain the bosonized 
expression of the interchain repulsion V±: 



V±a dx 



; \/2(0pj - ) cos \/20o-i cos V20cr2 + 



To obtain the bosonized expression of the spin exchange interaction, we need the spin density operators 



1 

7ra 
1 

7ra 



— — cos \f2Q„ p{x) cos \/20cr p(x) 

sy{x) 



sin V26a:p{x) COS V2(l)cr.p{x 
i\/2ipp^p — 2kFX 



2'Ka 

^i\/2(t)pp — 2kFX 



■ COS 



2 1 e 

Try z 



27ra 



27ra 
sin\/2(?!)o-,p(a;), 



nnV29c,p{x), 



.(17) 

(18) 
(19) 
(20) 



where d^^p = tt J^Ili, p{x')dx'. Since the terms coming from the 4:kp harmonics are less relevant than those produced 
by the 2kp harmonics, in Eqs. (|18|I - H2U|) we have altogether neglected the Akp harmonics, whose expression can be 
found in the Appendix IbI From the bosonized expression of the spin densities, we obtain the exchange interaction as: 



H, 



Ja_ 
" 27r2 

_J±_ 
■(27ra)2 



dxdx4'a,ldx4>a 



J± 



2{nay 



■ / dxcos \/2(0.,i -^^,2)cosV2(0,,i 



dx cos \/2(0p,i - 0P.2) [2 cos V2{e^s - Oa.2) + cos \/2((/)<,,i - ^^,2) - V2{(f>^s + 0<T,2)j ■ (21) 



The total interchain interaction Hamiltonian H± = -ffvj_ + ^ more conveniently written by introducing the 
new canonically conjugate fields^^ (j)^^± = {4>y,i ± (1)^,2)/^ and Iiu.± ~ {^u.i ± ^v,2)l^, as: 



+ 



(4V^ + Ji_)a 

[2-Kaf 
2CPV^ 
(27ra)2 



(27ra)2 



[ 2Jj_a f 

/ dx cos 20p_ cos 20cr+ + 7 TTT / dx cos 20„_ cos 29cr- 

J (2TTay J 



J dx cos 20p„ cos 200— 
J da;(cos40p_ + cos40p+) + J dx [{dx<j)p+f - {dx(l)p-)^] + ^ J [i^'- 



- {dx(l>c 



22) 



From the bosonized expression (|22|l of the interchain interactions and the Hamiltonians of the decoupled chains 
(|14|l - (|15ll we obtain the full bosonized form of the Hamiltonian (Ql-©- The total charge mode (j)p+ decouples from the 
spin and staggered charge modes. Below, we will start discussing the properties of the Hamiltonians of these modes. 



B. total charge Hamiltonian 



The total charge excitation Hamiltonian is 
where: 



Up+Kp+inUp+f + ^{dxcb, 



K. 



•p+) 



p+ 



Up[l + 



Kp+=Kp\l 



go = C^Vi_ 



2KpV±a 



2KpVj_a 



(23) 

(24) 

(25) 
(26) 
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The expression of the interchain couphng go in H23|l can also be derived in perturbation theory by the approach 
of Ref. [ij and a sketch of such derivation can be found in the Appendix ^ We note that the same sine-Gordon 
Hamiltonian was derived in Sec. Ill Al in the Umit U/t ^ oo, with the identification (p^ — \/2(f)p, K = 2Kp, — 2Kp^, 
and C = 1. A similar Hamiltonian was also obtained in Sec. Ill Bl in the limit J± oo. This fact indicates a continuity 
relation between the weak and the strong coupling regimes for the charge excitations. Finally, the Hamiltonian for 
the total charge mode H23|) can also be recovered from a general argumen1ii2i2^: From Eq. (|15|l . it is easily seen 
that a translation by one lattice site amounts to making \/2</'p,n — ^ \/20p,n — ""/S, and thus (f>p+ — > (f>p+ — 7r/2. 
Translation invariance of the lattice Hamiltonian requires the continuum bosonized Hamiltonian to be invariant under 
such transformation. The most relevant Umklapp operator compatible with this symmetry is cos4(/)p_|_, and it should 
be the operator responsible for the opening of the gap in the charge modes at quarter-filling. As a result, in the strong 
coupling regime, we find a boundary between the gapped and the gapless regime given by Kp^ — 1/2. Concerning the 
intrachain Umklapp terms, they are of the formlSiii cos 4-^/201,2 and become relevant only for Kp < 1/4. Thus we can 
safely neglect them compared with interchain Umklapp terms. A final remark is in order concerning the sign of go. In 
appendixiBl we show that the interaction in the spin sector only produces a term of the form cos 4:(/)p+ cos 20„- cos 20^+. 
This latter term is less relevant than a term of the form cos4(^p_|_, so that we can usually neglect it and assume 
go > 0. However, if the spin gap in the system is large, this term cannot be neglected and contributes a correction 
cos4(/)p_|_(cos20cr- cos2</)cr+) to the Hamiltonian (|23|) . This can result in a change of the sign of go. In that case, 
instead of having {<f>p+) = |- in the ground state, we should have instead ((/)p+) = 0. In the following we will discuss 
both the phases with (0p+) = j and {(l)p+) = 0. The latter ones are expected to be obtained when J± ^ V± since we 
need a large spin gap to modify the sign of go. 



C. spin and charge difference Hamiltonian 



From Eqs. H14|) - (|15|l and Eq. H22|) . we obtain the Hamiltonian describing the interaction of spin modes 0o-± and 
interchain charge modes (/)p„ . Under renormalization group (RG) transformation, besides the interactions already 
present in the bare Hamiltonians, new interactions can be generated. Thus the Hamiltonian to consider reads: 



H = 



where: 



251 



ue{p—,(T+,cr—} 



(27ra) 



,„ ,„ , da; cos 2(/)cr+ cos 20cr- i 
(27ra)^ J (zTTa)'^ 

/ dx cos 20p_ cos 20^+ -I- — — — 



232 



dx cos 2(j)cr+ COS 26a- 



2,93 



(2™) 



dx cos ' 



/ dx cos 20p_ cos 20£r- + T— / da;cos2(/)p_ cos26'ct-, 
J (2iTay J 



2Vj_Kpa 



1/2 



Ua+ = Ua 



2TrUn 



1/2 



1 - 



2nvF J 



1/2 



Kn 1 



21/_L ifpfl 



-1/2 



-1/2 



9i = 0, 52 J^a, 53 = 2C^V^a, 



Ji a 



-1/2 



2ttvf I 



J±a. 



(27) 

(28) 

(29) 

(30) 

(31) 
(32) 
(33) 



Note that in Eqs. H28|) . we have put gi ^ because the intrachain spin interaction in the repulsive Hubbard model are 
marginally irrelevant. In full rigor, we should have added a small marginally irrelevant interaction to the Hamiltonian 
(|15|l and we would have obtained 51 > and slightly modified expressions of Kcr± . In spite of the fact that we have 
put 51 = in the initial conditions, in the following we will also discuss the case of a relevant 51 under RG. 

To study the Hamiltonian H27|) RG equations can be derived. Velocity differences between the various modes can 
be neglected, and one can take Ua-± — Up- = vp- Introducing yi = gi/{-KVF) and K„r = 1 + y<jrl2, and using the 
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operator product expansion (OPE) niethods2&, the RG equations read: 

= yl + \{yl + yl + yll (34) 

\{yl + yl-yl-yl), (35) 
1 



-( 


1 


dl \ 


Kp- 


dya- 




dl 




dya- 


f _ 


dl 




dyi 




dl 




dy2 




dl 




dy3 
dl 


-(2 


dyi 
dl 


= (1 


dys 
dl 


= (1 


dye 


= (1 


dl 



' ' yl + yl), (36) 

[(2/<T+ + ya-)yi + y^ys] , (37) 

[(y<T+ - ya-)y2 + y^ye] , (38) 

4i^p_)y3-i(y|+2/| + %'), (39) 

Kp-)y4 ~ ^{ya+y4 + yiyz + y2ye + 2/3^4), (40) 

Kp-)y5 - 7:{ya~yb + ysys + yiyi), (4i) 



(42) 



We see from these equations that although yi is initially zero, it becomes non-zero under the RG flow. The fact that 
the problem under consideration has SU{2) spin rotational invariance leads to a simplification of the RG equations 
(|34ll . The initial conditions H28|) lead to the following relations Vl: 

y,+ {l)+y,^{l)=y,{l), (43) 

ya-{i)-y<r+{i) = y2{i), (44) 

2/5(0-2/4(0=2/6(0- (45) 

These conditions ensure the SU (2) symmetry of the RG flow and reduce the flow from a curve in a nine-dimensional 
space to a curve in a six-dimensional hyperplane. The simplified RG equations now reads: 



^ = -l{y! + yl + y4ye) (47) 

^4(2/.^-2/42/a) (48) 

^ = (2 - 4Kp^)y, - \{yl + yl + y^ye) (49) 

^ = (1 - Kp-)yi - ^(32/12/4 " 2/22/4 + 22/22/6 + 22/i2/6 + 22/32/4) (50) 
dye / \ 1 / 

-rr = (1 - Kp-)ye + 7(2/12/6 + 2/22/6 - 22/32/6 - 22/22/4)- (51) 



The flow of these equations will be discussed in Sec. IVII Here, we want to discuss briefly the possible fixed points 
within a semiclassical argument, i.e. by considering the expectation values of the phase fields that minimize the 
classical ground state energy. By looking at the possible classical minima, we can distinguish two regimes. In the 
first regime, {4>p-) =0,^ and (cos20p_) ^ 0. This case is similar to that obtained in ladders at incommensurate 
£]^jjj^g30j3ij74 ^jjjg regime the relevance of the terms (74,5,6 is responsible for the presence of the spin gap. In the 
second regime, {(t)p-) = so that (cos40p_) ^ 0, and (cos2(/)p_) = 0. This regime corresponds to 33 being the 
dominant interaction, and (74,5,6 being less relevant interactions. This case is specific of the quarter-filled ladder, and 
corresponds in the limit ?7 = c» of Sec. Ill Al to the formation of the gap in (/)_. We note that in the present regime, 
spin gaps can also be generated by the (71,2 terms. When gi is the relevant interaction, the corresponding bosonized 
expression can be rewritten as ~ gi (cos •\/8(/)i -I- cos ^/S(f)2 ) , so that the spin gaps correspond to two independent 
intrachain spin gaps. We also note that in order to preserve spin rotational invariance, we need to have gi < in that 
case. 
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It is possible to give a roughly estimate of the parameters region where each regime dominate by looking at the 
scaling dimensions of the operators. The operator cos4(/)p_ has scaling dimension 4:Kp_ so that it is relevant for 
Kp^ < 1/2. The operators cos20p_ cos20o._, cos2(/)p_ cos2(/)o-+, cos2(/)p_ cos 26*0— have respective scaling dimensions 
Kp- + K„-, Kp- + K„-\., Kp- + 1/ Kc-. Taking into account the spin rotational symmetry, this implies that these 
operators have all the dimension Kp^ + 1 and become relevant for Kp^ < 1. The operator cos4(/)p_ becomes the 
most relevant operator when 4Kp- < 1 + Kp-, i. c. for Kp- < 1/3. Therefore, we expect to have the first regime in 
the limit of a repulsion not too strong, 1/3 < ifp_ < 1/2, and the second regime in the case of a stronger repulsion 
1/4 < Kp_ < 1/2. 



IV. PHASE DIAGRAM 



In this section, we describe the various insulating phases predicted from the renormalization group equations H51|l . 
In order to describe the possible ground states in the phase diagram, we introduce first the corresponding order 
parameters. 



A. order parameters 



Since for strong repulsion the system has a gap in the spin excitations, the possible order parameters can 
only be bond-order waves (BOW) and charge density waves (CDW). We will denote the corresponding phases as 
{QxiQy)~^OW, and (g^,, g^)— CDW where = tt/2 or tt and (/j, = 0, tt or Qy = ±7r/2 if = 7r/2. The phases with 
a qx — ordering have at least two-fold degenerate ground state, while the phases with a qx = ^ have a four-fold 
degenerate ground state. 

On the lattice, these operators are defined as: 

OcDW(,T,qy){-i) = + cos{qy)n^^2){qy = 0,7r) (52) 

cr 

OcDW(f ,,,)(*) = e-'^^(e*«»/2^.,i +e-*«»/2n,,2)(g, = 0,|,7r,-|). (54) 
The corresponding bosonized expressions is obtained by those of the charge densities as a function of (j)p± and (t>ai,2- 



1 „H<Pp++<Pp_-2fcF2:) p 

Plix) = — a.(</>p+ +0pj + cosV2(f>,,{x) + }:_eM^.^+4>,_-2k,x)^ ^^^^ 

TT 7ra Tra 

p,{x) = —dx{(Pp^ - <t>p^) + - cos V20,,(x) + ^^^^i^.+ -4'.-~2k,x) ^ (5g) 

TT Txa na 

We consider first q^ — n. The charge density wave order parameters are straightforwardly obtained from (|55|I - H56(I as: 

OcDWiTTfi) ^ cos 20P+ cos 20p_ (57) 

OcDWiTT^TT) ~ sin 2(?!)p+ sin 2(/)p_ (58) 

The bond order wave order parameters measure the charge density between the sites i and i + 1, i.e. on site i + 1/2, 
and their bosonized expression reads: 

Obow(,t,o) ^ sin 20P+ cos 2(/)p^ (59) 

OBOWirr^TT) ~ COS 2(?!)p+ sin 2(/)p_ (60) 

Concerning the tt/2 charge density wave order parameter, we have to consider whether (j)^- or is ordered. In the 
case ^CT- is ordered, the operators cos •\/20o-,i,2 have zero expectation value, and exponentially decaying correlations, 
so that all of the {tt/2, qy)-CDW order parameters have short-range order. When ^o-- and (pa+ are ordered, the 
operators cos \/20(t,i.2 have non-zero expectation values. In this case, the expression of the order parameters for 
{tt/2, qy)— CDW as a function of 4>p+ and 0p_ is: 
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(0P-> 




(<^-> 




(f ,7r)-CDW 
(7r,0)-CDW 


±f 








- 


(7r,0)-CDW 


±f 










(7r,7r)-BOW 


±f 










(f,f)-CDW 
(7r,7r)-BOW 


±f 


2 








TABLE I: The phases of the quarter-filled ladder in the case (<^p+) = 





(0P-> 








(f ,7r)-CDW 
(7r,0)-BOW 


±f 










(7r,0)-BOW 


±f 










(7r,7r)-CDW 


±f 










(f,f)-CDW 
(7r,7r)-CDW 


±f 











TABLE IL The phases of the quarter-filled ladder in the case {4>p+) = -f 



OcDWi^^O) - ^6^"^"+ COS (/.p. (61) 

OcDWi^.^) - ^e^'^"^o«(<^p--J) (62) 
OcDWi^,.) - ^e^-^^+sin^p- (63) 

OcDWi^^-i) - ^^''''^<^os{c^p- + l) (64) 

From the knowledge of the order parameters, in the following we are going to discuss the ground states that are 
realized for the various possible ordering of the phase fields. We will first discuss the orderings with {(j)p-) = j 
and then turn to orderings with {(j>p-) = 0. As we have explained in Sec. IIII Bl the latter type of ordering should 
be expected for dominant J±. In Tables llllll we give a summary of the phases corresponding to {4>p+) = 7r/4 and 
{(t>p+) ~ 0. These different phases are detailed in the forthcoming sections. 



B. (tt, 7r)-charge density wave 

In the case of (0p+) = {4>p-) = j, {4>(y+) = 0, {6^-) = f, the (7r/2, gj,)— CDW order parameters have all 
zero expectation value and exponentially decaying correlations. Moreover, the (tt, 7r)-CDW order parameter has a 
non-zero expectation value, leading to a zig-zag charge ordering^i*^. This ordering is represented on Fig. ^ This 
charge ordered phase possesses also a spin gap caused by interchain coupling and identifies with the CDWgg (with 
spin gap) of Ref. 48. The ground state is this phase is fourfold degenerate (twofold due to translation symmetry, 
and twofold due to the two possible orientations of the spin singlets). We note that this phase has previously bee 
discussed in Sec. Ill Al in the limit oi U — > oo, however in that limit it was not possible to discuss the spin modes. 
From the discussion of Sec. Ill Al we see that the essential ingredient for the zig-zag ordering is the mutual locking of 
the 4:kp charge density fluctuations so that zig-zag ordering needs both strong intrachain and interchain repulsion. 
The formation of the spin gap appears to be unrelated to the zig-zag ordering, but only a consequence of the coupling 
of the q ~ fluctuations. We note however that the higher order terms derived in App. could lead in the case 
of a large charge gap to corrections to the spin Hamiltonian that could enhance the spin gap. In the RG study of 
Eqs. (|34|l . this phase is obtained for j/3 — > +00, and y2 — *■ -l-oo. 
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FIG. 1; The phase with zig-zag charge ordering and spin gap- (7r,7r)-CDW. The grey circles represent the local charge density 
and the arrows the spins. The empty circles indicate the absence of spin or charge on that site. 

C. (f if) Charge density wave 

In the case of (0p+) — {4>p-) — ±f and {(})„+) = 0, {<pa~) = a (|,±f )-CDW state is formed coexisting with a 
(tt, 7r)-CDW oscillation. This phase possesses also an intrachain spin gap. The corresponding state is represented on 
Fig. 13 A simple physical picture of this state is that a ^--charge density is formed in each chain, and the interchain 
repulsion locks the phases of both charge density waves. Since there is an intrachain gap, this phase appears to be 
more likely to be observed in a Hubbard ladder with [/ < and V|| > or in a t-J ladder in a regime in which Jy 
is large enough to cause the formation of a spin gap in the single chain*"^. The dephasing between the two charge 
density waves is 2(0p_) = ±f . This dephasing results from the mutual locking of the Akp density fluctuations, thus 
pointing to strong repulsion in the chain. The corresponding phase has a fourfold degenerate ground state. In the 
rcnormalization group treatment, this phase is obtained for yi — > ~\-oo and 2/3 — > -l-oo. 

FIG. 2: The (-1, -D-charge density wave phase. The grey circles represent the local charge density. The difference in size 
indicates small or large charge density on that site. 

D. (tt, 0) Bond order wave 

For (0p+) = f , {4>p-) = f , {4'(y+) = Oj (^CT-) = 0, all of the CDW order parameters vanish. The only nonvanishing 
order parameter is (tt, 0)-BOW. This BOW can be described in physical terms in the following way: the fermions are 
localized on the bonds between two sites, and due to the Jj_ interaction, they form a spin singlet with the fermion 
on the opposite chain. Clearly, this phase should be expected at strong Jj_ and moderate repulsion in the chains. 
The ground state is only twofold degenerate. In the rcnormalization group, this phase is obtained for yi —t —00, 
2/5 + 2/6 ^ +00, and the inspection of the expression of y^ confirms that Jl is the dominant interaction in the (tt, 0)- 
BOW. The corresponding phase is drawn on Figure |21 We note that since this phase has a uniform charge density, it 
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FIG. 3: The (tt, 0)-bond-order wave phase. The connected black circles represent a bond occupied by an electron. 



is a possible candidate for the homogeneous insulator Hlgg phase of Ref. 48. It is also interesting to remark that this 
phase is a two-leg analog of the SDWi phase obtained irJI- in a system of coupled quarter filled extended Hubbard 
chains. 

E. Charge density wave 

For {(j)p+) — J, {<j)p~-) — f , (0(7+) = 0, = 0, the (5,7r)-CDW order parameter does not vanish. The (7r,0)- 

BOW correlations are also present in this phase. A sketch of this phase is given in Fig. 01 Similarly to the ■|)-CDW 
this phase results from the mutual locking of the 2kp charge density wave fluctuations of the two coupled chains. 
However, in contrast to the ^)-CDW this locking is produced by the coupling of the 2kF density fluctuations. 
This implies that the spin gap is formed as a result of interchain coupling. In the RG study this phase is obtained for 

M H ^ 




FIG. 4: The (-1, 7r)-CDW with (7r,0)-BOW. The grey circles represent the charge density in Fig|5|and the black horizontal 
lines denote occupied bonds as in Fig. |21 In contrast with the pure (vr, 0)-BOW the charge density is inhomogeneous along the 
chains. 



2/4 +00, 2/5 + ye ~* +00. Using the expression of j/4 we see that this phase is dominated by interchain repulsion in 
contrast with the (tt, 0)-BOW. 
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(tt, tt) Bond order wave 



When 



= 0, 



= ±- 



= 0, all of the CDW order parameters vanish and only the (tt, tt) 



BOW order parameter is non-zero. This phase corresponds to a staggered bond ordering shown on Fig. |S1 It can be 
viewed as a (tt, 7r)-CDW translated by half a lattice spacing. Elementary excitations in this phase are thus similar to 
those of the (tt, 7r)-CDW phase. In the renormalization group such state is obtained for y2 +oo, +oo. 




FIG. 5: Sketch of the (tt, 7r)-BOW. 



G. 



Charge Density Wave 



When 



±'- 



0, the (J,?)-CDW order parameter does not vanish and 



coexists with a (tt, 7r)-B0W. This phase is different from the (^,^)-CDW previously encountered as the previous 
phase contained a (tt, 7r)-CDW. There are however similarities between these two phases since both result from the 
mutual locking of Akp components of density fluctuations combined with an intrachain spin gap. Since the sign change 
in (7o results from the formation of an interchain spin gap, we should expect this phase to have only a rather narrow 
domain of existence. This phase is sketched on Fig. El In the renormalization group treatment, it is obtained when 
yi +00, ys -> +CX). 



4-t 



4 




FIG. 6: The (f , f )-CDW with (7r,7r)-BOW. The grey circles have the same meaning as in FigH 



H. (tt, 0) Charge Density Wave 



When (0,+) = O,(0p_) = ±f , (0,+) = f , ifi^J) 
one. In the RG study this phase is obtained for 2/4 



0, the only non-vanishing order parameter is the (tt, 0)-CDW 
—00, j/5 -l- ye +00, i.e. it corresponds to a dominant Jj^, in 
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agreement with the results of Sec. Ill Bl This phase can be viewed as the (tt, 0)-BOW shifted by a half lattice spacing 
The corresponding phase is sketched in Fig. [7| 



O 



O 




O 



O 




FIG. 7: The (tt, 0)-CDW. The grey (empty) circles have the same meaning as in FigQ 



I. (f ,1") Charge Density Wave 

When ((^p+) — O,(0p_) — ±-|, (0cr+) = f , {<t>a-) — 0, the only non-vanishing order parameters are the (7r,0)-CDW 
and the (|, 7r)-CDW ones. Under RG study, this phase is obtained for y4 +oo, 2/5 + ye ^ +00 and corresponds to 
dominant VI. For this reason, we should expect this phase to have a rather narrow domain of existence. The spin 
gap corresponds again to an intrachain spin gap. This phase is sketched in Fig. |S1 






FIG. 8: The (|,7r)-CDW with (tt, 0)-CDW. The grey circles have the same meaning as in FiglH 



ELEMENTARY EXCITATIONS 



In the present section, we discuss the nature of the elementary excitations in the various gapped phases of the 
quarter-filled two leg ladder. Due to the spin-charge separation, all of the insulating phases have charge excitations 
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of charge ±e and zero spin. However, the nature of the magnetic excitations is dependent of the nature of the phases 
considered. Below, we review the various phases with their elementary excitations. 



In these two phases, - 
sine- Gordon Hamiltonians: 



A. (tt, 7r)-CDW and (tt, 0)-BOW 
1. Charge excitations 

-oo, implying that the effective Hamiltonian for the charge sector is described by two 



where 



Hcff - Hi 



-{dx4>p 



/ rf.cos.,,„ 



(65) 



(66) 



a = ±. Thus, the elementary excitations can be described in terms of solitons of two decoupled sine-Gordon models. 
Following a semiclassical argument, a soliton is joining two consecutive minima of the Hamiltonian 1)66(1 . These 
minimums are given by < (jyp^ >= 7r/4[7r/2]. As a result, the charge of solitons is given by: 



+ 00 



a 



[(/.pJ+«3) - 0pJ-C5o)] = ±1 

TT 



(67) 



These solitons can be understood as domain walls between the two charge ordered ground states as represented 
on figure IHl^a) and figure |5{c). Alternatively, the charged solitons can be viewed as holon-holon (for charge — e) or 
antiholon-antiholon (for charge -l-e) bound states. The neutral solitons can similarly be viewed as holon-antiholon 
bound state. An immediate consequence of the mapping of the charge excitations onto the sine-Gordon model is 
that the transport properties of the quarter-filled ladder can be obtained from the methods reviewed in Ref.®^. For 
Kpa > 1/4, the solitons are the only possible excitations in the model. For more strongly repulsive excitations, the 
formation of soliton bound states (breathers) becomes possible. The existence of these excitations translates into an 
exciton peak in the optical conductivity. Such exciton peak has been studied in the case of the single Hubbard chain 
in Ref. |73- However, in the model we are considering, for Kpa < 1/4, intrachain Umklapp scattering becomes relevant 
and total/antisymmetric charge excitations no more decouple in this regime. 



2. spin excitations 



To describe spin excitations, we need to consider the effective Hamiltonian of the spin modes: 

2.g2 



r=± 



-idx4>o 



(27ra) 



dx cos 2(/)ct+ cos 26a 



(68) 



This Hamiltonian corresponds exactly to the bosonized Hamiltonian of a zig-zag spin ladderiSiffiiSi. Remarkably, in 
the strong coupling limit VpV±,U t±,t, a mapping on a zig-zag spin chain was derived to describe the low energy 
excitations in the (tt, 7r)-CDW phased. We thus notice the continuity between the weak coupling and the strong 
coupling limit for the spin excitations of the (tt, 7r)-CDW in this problem. Spin excitations of the quarter filled ladder 
are those of a zig-zag ladder, i. e. massive spinous. A more detailed picture of the spin excitations can be gained 
by applying a transformation due to Witten and Shanka^SiiSSiSi. We first notice that we have Ka+K^- = 1 and 
Ua-+ = Ua-. Thus, we can perform a duality transformation (j)^- = 9a-,0a- = ^ct-, and introduce: 
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FIG. 9: Excitations of the quarter filled ladder in the zig-zag charge ordered phase: 

(a) Charged excitation: it can be viewed as a domain wall between the two zig-zag charge ordered ground states of the ladder 
or equivalently as a bound state of two antiholons of the quarter-filled chains. 

(b) spinon excitation. It can be viewed as a domain wall on the spontaneously dimerized effective zig-zag ladder or as spinon 
inside one of the chains. 

(c) Neutral excitations. It can be viewed either as a domain wall between the two zig-zag charge ordered ground states or as a 
holon-antiholon bound state. 



This procedure permits us to reduce the effective Hamiltonian to: 



"'-It. 



UbKb{T:Ilb 



92 



(27ra)2 

92 



dx cos 
dx cos 



(71) 
(72) 

(73) 



where Ua = = Ua+, Ka = Kb = Ka-+ and 52 preserves SU{2) symmetry. The original Hamiltonian Ht)8|l is thus 
decoupled into two massive sine-Gordon model at the SU{2) point. They describe two massive spin excitations 
carrying spin 1/2. In the strong coupling limit, the spin excitations can be viewed as the spinous of a zig-zag ladder 
(see figure 



B. (f,f)-CDW 



In the (^, ±^)-charge density wave phases, the charge excitations carry the same quantum numbers as in the zig- 
zag charge ordered phase. They correspond to domain walls between the four different (^, ^)-CDW ground states. 
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However, spin excitations are of a different nature since both Sf and S'| are good quantum numbers. In fact, these 
spin excitations are massive spinous "confined" within each chain. 



C. phases with {4>p-) = 

The total charge excitations still carry ±e, as in the previous phase, but this time it is not possible to decouple spin 
excitations from antisymmetric charge excitations. When (p^- is ordered, we have (/)p_(+oo) — (j)p-{—oo) = ±7r/2, 
(/)o-+(+oo) ~ '/'ct+(^oo) = ±7r/2 and 0(j-(+oo) — (— oo) = ±7r/2. As a result, the elementary excitations carry total 
charge zero, charge difference ±e, and spin SI = ±1/2 or 5| = ±1/2. These elementary excitations can be viewed 
as a three body bound state of a holon in one chain, a antiholon in the other chain and a spinon. This excitation is 
sketched on Fig. ^1 



A 



A 




(a) 



W 




A 
Y 

A 



V7 



(b) 



FIG. 10: Elementary excitations in the (tt, 0)-CDW. (a) holon and antiholon excitations carrying the charge ±e respectively, 
(b) Mixed spin/orbital excitation. 



In fact, there exists an interesting analogy between the spin/antisymmetric charge sector of the quarter filled ladder 
and the half-filled Hubbard-Kondo-Heisenberg (HKH) chain. 
The half-filled HKH chain is described by the Hamiltonian: 



H 



-t 



„t 



2 



E 



(74) 



where a^'^'^ are the usual Pauli matrices. At half-filling, the resulting bosonizcd Hamiltoniai 



.84.85.86 



reads: 



i7 = 



dx 
2^ 
dx 
2^ 



dx 
2^ 



2Jk 



{dx(l>o 



(27ra)2 



dx cos a/20p [cos 20O-- 



(75) 



cos 20^+ + 2 cos 26'^_] (76) 



Making a rescaling 0p = ^/2(f>p-,Kp = 2Kp-, we obtain the same Hamiltonian as the one describing the spin/charge 
difference excitations of the quarter-filled ladder H22() - (|27|l . In particular, the point Kp- — 1/2 in the quarter filled 
ladder corresponds to Kp = 1 in the Kondo-Heisenberg chain, i.e. C/ = in (|74|) . Thus, magnetic properties near the 
metal insulator transition of the quarter filled ladder should be analogous to those of a Kondo-Heisenberg insulator. 
Close to Kp^ — 1, the structure of excitations is quite different. In the vicinity of that poinli^I, gapful excitations of 
the ladder away from half-filling are described by a massive 5*0(3) x 5*0(3) Gross-Neveu modcl'*^ as a consequence 
of the breaking of the original 5C/(4) symmetry to SU{2) x 51/(2). The enhancement of the symmetry by the RG is 
instead absent for A'p_ = 1/2. 
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VI. NUMERICAL STUDY OF THE PHASE DIAGRAM 

In this Section we discuss the numerical results of the RG equations. Prior than that we analyze in details the nature 
of the transition among the various phases. The analysis is partially accomplished by a mapping onto a Majorana 
fermion theory, that gives also rise to a connection with the spin-orbital models. 

A. Mapping on a theory containing Majorana Fermions 

We consider first the part of the Hamiltonian H27() with couplings 54. 5.6- It is convenient to rewrite this part of the 
interaction in terms of massive Majorana fermion operators {v = 0, 1, 2, 3). Using the identities'^: 

coa2<j)^+ 

= i[Cr,iCl.1 + Cfl,2CL,2), (77) 



na 
cos 2(/)ct 



j(Ci?.3a,3 + Ci?.oa,o), (78) 



7ra 

and the relation 35 = 34 + ge: the interaction terms 54,5,6 are rewritten as: 



na 

cos20^_ /> > > > X 

= «(Ci?,oa,o - Cfl.aCL.a), (79) 



i cos 2(/)p_ 



7r^(Cfl,iCi,i + Cfl,2Ci,2 + C-R.sCi.a) + ^ — —Cr.oCl.o 

2na Zna 



(80) 



cos^ 


9<T+ = 0'icr2, 


cose; 




sinq 




cos 


da- = AisCTo 



where 175 + is always a positive quantity. 

In the continuum limit, following ReffSSI we can express the bosonic exponents in terms of the order (ai) and 
disorder (/i^) parameters of four Ising models as follows: 

(81) 
(82) 
(83) 
(84) 

When {4>p-) = ^, the Ising order parameter cto has a non-zero expectation value. However the sign of 34 is not 
fixed. For 54 > 0, and {4>p-) = ^, ci,2,3 have all nonzero expectation values, corresponding to have both and cj)^- 
ordered. For 54 < 0, /ii.2,3 have all nonzero expectation values, corresponding to having both ipa-^ and 9a-- ordered. 

Now, we would like to focus on the interaction part of the Hamiltonian that contains the terms ffi, (72,ffo-±- Using 
the mapping on Majorana fermions, we can write this part as: 

^(9,0,+ )^ + ^(a,0,+ )2 + -^I^cos2(/),_cos20,+ + --/i_cos20,_cos20,+ (85) 
27r^ 27r^ 2(7ra)^ 2[na)'^ 

= ^9(T+C-R.,lCi,lCfl,2C-L,2 — .g(T-Cfl,oCL,oCi?,3CL,3 " 5+(C-R.,lCi,l + C-R,2Ci,2)Cfl,oC-L,0 " .9- (C-R,lC-L,l + Cr,2Cl,2)Cr,3Cl,3- 

where g± = {gi ± g2)/2. Using the SU{2) symmetry conditions l|43(l . this expression can be further rewritten as: 

~ 5cr+(C-R,lCL,lC-R,2CL,2 + Ck,iCl,iC-R,3Cl,3 + C-R,2Ci,2Cfl,3CL,3) - .9(t- (Cfl,lCL,l + C-R,2Cl,2 + Cfl,3CL,3)C-R,oCL,0, (86) 

which makes the SU(2) symmetry transparent. We note that for g^^ > 0, i(Cfl.aCi,a) can have opposite signs 
depending whether a = or a 7^ 0. This implies that either (^0) 7^ and (0-1,2,3) 7^ or (cro) 7^ and (^^1,2,3) 7^ 0. 
Both cases correspond to having 9^^ and (f)a--i- ordered. When g^- < 0, we have both 0^+ and (j)^- ordered. 

1. Spin-orbital models 

Before closing this section, we would like to discuss the following connection. If we consider the case of V\\ = 0, 
J± ~ and V± ~ U, we can rewrite the interchain interaction in Eq. |^ as: 

^ + n,A,i + nj,2,T + "-^,2,i)^ (87) 
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which has a manifest 5/7(4) invariance. Thus the problem is related to the quarter-filled 5't/(4) Hubbard model^S. 
The charge Umklapp term derived for that modelSS agrees with (|23f) . If U is large enough, the low energy excitations 
are those of an antiferromagnetic chain of SU{A) spinsSfi. Considering deviations from the S'C/(4) symmetric point, 
it is possible to derive a model describing the low energy excitations of the insulator^^ in terms of coupled spin 
and orbital modes (the spin-orbital model). In Refs. ,.88i.91...92 . a SU{2) x SU{2) spin-orbital symmetric model was 
analyzed perturbatively around the SU{A) point using bosonization and rcfcrmionization techniques. An Hamiltonian 
describing the low energy dynamics of the system with six Majorana fermions was obtained, and the formation of a 
dimerized spin gapped phase was predicted in the physical range of parameters. The existence of this dimerized phase 
was also confirmed by numerical studies22i2^. In our present problem, assuming Kp^ ^ 1 so that the term containing 
f/3 can be neglected, and using the mapping: 

^^^^ = ^{Cb.aCla + Cb.aCla), (88) 
na 

we can recast the interactions in the spin/antisymmetric charge sector of Eq. (|27|l in terms of two triplets of Majorana 
fermions: (Ci, C2, Cs) representing the spin excitations, and (Co, C4; Cs) representing the orbital excitations as in R.ef . Is^ 
The Hamiltonian we obtain is however more general, since we did not assume any SU (2) symmetry in the 
orbital sector. From our previous discussion of the phase diagram in Sec. IIVI we expect that the dimerized phase 
(Phase IV in Ref. or Phase A in Ref. IqH) of the spin-orbital model is related to the (tt, 0)-BOW. In a different 
limitSMi of the SU{2) x SU{2) spin-orbital model, a dimerized insulating phase analogous to the (7r,0)-BOW was 
also obtained. We note that no charge ordering is predicted by the SU{2) x SU{2) spin orbital model. The reason 
for this is that the V|| term of the ladder model produces a large renormalization of Kp^ in the charge ordered state, 
while leaving K„-, Ka-+ close to the non-interacting value, due to the remaining SU{2) symmetry. To describe charge 
ordering, one has thus to consider more general spin-orbital models in which the interaction in the orbital sector has 
only the U{1) symmetry^. 



B. order of the transitions between the different phases 



We now turn to the quantum phase transitions between the different phases of the quarter-filled ladder. These 
phase transitions can result from a change of the ordering in the <pp- field, the spin gap being preserved, or from a 
change of the ordering in the spin sector. As we will see below, the former phase transitions are of second order and 
belong to the one-dimensional quantum Ising universality class, whereas the latter transitions can be either of second 
and first order. We begin discussing the Ising transitions in the antisymmetric charge mode. 



1. Ising transitions 



We consider phase transitions in which the order in the spin sector is not changed, but the order in the antisymmetric 
charge sector (/>p_ is modified. These transitions occur between the (tt, 0)-BOW and the (tt, 7r)-CDW, the 7r)-CDWs 
and the (^, ^)-CDWs, the (tt, 7r)-B0W and the (7r,0)-CDW. Since the order in the spin sector does not change we 
can describe the transition by concentrating only on the (j)p^ modes. The resulting effective Hamiltonian reads: 



iip_ii'p_(7rHp_) 



233 



(27ra)2 



dx cos ' 



25 



(27ra)2 



dx cos2(/)p_. 



(89) 



This Hamiltonian is the one of the double sine-Gordon model^®. The semiclassical analysis of the double cosine 
potential shows that for ((/ip-) = ±j, (cos2(/)p_) = 0. Using the result of^-^, this implies that as g is varied, an Ising 
transition is obtained in 4>p-. For the sake of simplicity, let us consider the transition between the (7r,0)-BOW and 
the (tt, tt)— CDW. Similar features appear in all the other Ising transitions. In the (7r,7r)-CDW, the order parameter 
OcDW(„ ^ sin20p_. Using the results o^ on the ultraviolet-infrared transmutation of operators, it is easily 
seen that OcDW^■„ „) ^ where ^ is the disorder parameter of a quantum Ising model. This implies in particular 

that {OcDW^„ ^){x)OcDW^, x) (0)) ^ x~^^'^ at the transition. Moreover, near the transition one has: {OcDW^„ „)) ^ 
1/8 

. If the gap Ap_ vanishes linearly with the interaction V this imphes (Oc_dW(„.„)) ~ {V —VcY'^ , giving rise 

to the onset of the order parameter at a critical value Vc- Such onset was numerically observed at the homogeneous 
insulator HIjg-CDWsg transition in Refli^ and attributed to an Ising transition. In our results, we observe that 
neither the charge gap nor the spin gap go to zero at the transition, meaning the absence of a critical value Vc for the 
onset of the transition. 
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A simple picture of the transition is obtained, in the limit of strong coupling Vj_ oo, by an effective quantum 
Ising model. In this picture, the only states retained are the electron pairs forming singlets along the diagonal and 
pointing either in the northwest or northeast direction. The variable — 1 when the diagonal singlet formed of two 
electrons is oriented northwest (NW), and cr^ = — 1 when the singlet is oriented northeast (NE). In the ground state 
all the electron pairs must have the same orientation. The corresponding potential energy reads: 

Fpot. =-^|^«+i- (90) 

n 

The kinetic energy comes from the term t. In second order perturbation theory, t flips a singlet pair from the NW to 
NE orientation. This process leads to a kinetic term in the Hamiltonian: 

^2 

Hkin. = (91) 

When the kinetic term in the Hamiltonian dominates, the singlet pairs go back and forth between NE and NW 
orientation leading to a zero average of the (tt, 7r)-CDW order parameter, and an effective bond order wave. When 
the potential energy dominates, the singlet pairs are all locked in the NE or NW orientation giving rise to a nonzero 
(tt, 7r)-CDW. 

2. spin transitions 

In the present section, we discuss the transitions in the spin sector. Since the gaps in the0p_ and (/)p_|_ sectors 
are robust, we can concentrate on a low energy effective spin model. First, let us focus on the case of (0p-) = §. 
In that case, the discussion is equivalent to the one in Ref. l52l The theory describing the transition point is the 
0(3) Gross-Neveu model, and the operator that takes the system away from the transition point is the mass of the 
Gross-Neveu fermions. As a result, the system will have a second order phase transition in the SU {2)2 WZW model 
universality class when the 0(3) Gross-Neveu has no spontaneous symmetry breaking and a first order transition 
when the Gross-Neveu model presents dynamical mass generation2£'2£. For ga+ < 0, there is a spontaneous symmetry 
breaking and thus a first order transition. Since this corresponds to J± > 0, first order spin transitions should be 
generic in the models we are considering. In particular, first order transitions should be expected between phases such 
as the (7r,0)-BOW and the (|,7r)-CDW or the (tt, 0)-CDW and the (|,7r)-CDW. These first order transitions occur 
in the spin sector and should be observable by looking at spin-spin correlations. In the case of {(j^p-) — j, we have 
to focus on the terms coming from 51 , 172 , ga± ■ In that case, for a transition to be possible we must have g^- = as 
one can see directly from Eq. H85|l . The behavior at the transition then depends on the sign of gcr+- When ga+ > 0, 
no gap is generated in the triplet modes thus giving a SU{2)i x SU{2)i criticality. For grr+ < 0, the triplet modes 
remain massive at the transition, leading to an Ising criticality. Since ga+ ^ — J_l < 0, an Ising transition is obtained 
between the (7r,7r)-BOW and the (f , f )-CDW or the (7r,7r)-CDW and the (f , f )-CDW. 

C. RG calculation 

To find the phase diagram, we integrate the Eqs. H51|l numerically using a fourth order Runge-Kutta algorithm for 
fixed values of Kp_ at varying Vj^ and J±. We stop the numerical integration when one of the coupling constants 
VsjUb^Ue becomes of order 1. We find that at this scale, j/1,1/2 are still inferior to 1. We have the following results. 
First, for A'p_ larger than 1/3, {Kp^ ~ 1/2) the ■|)-CDW and (tt, 7r)-CDW are absent. This is a consequence of 
the fact that in this regime, the 2fci? fluctuations are dominant over the Akp ones. As a result, we find either the 
bond order wave (tt, 0)-BOW or the (|-, 7r)-CDW. As could be expected, a large J± favors the former, and a large V± 
favors the latter. A first order transition is expected between these two phases. The phase diagram for Kp- = 0.5 is 
drawn on Fig. ^] When Kp- becomes smaller but still larger than 1/3, the term cos40p_ is more relevant and the 
(■|, |-)-CDW phase becomes stable at large V±. This can be understood as resulting from an increase in the strength 
of Akp fluctuations. The phase diagram for Kp = 0.35 is sketched on Fig. 1121 For Kp^ close to 1/3, the topology of 
the phase diagram becomes more complex. The (7r,7r)-CDW phase is present for competing V±,Jj_, along with the 
(tt, 0)-BOW and the two other CDWs. The phase diagram is sketched on Fig. ^1 Finally, for Kp^ < 1/3, cos4(/)p_ 
is the most relevant operator, and the (^,7r)-CDW phase disappears altogether. The (tt, 7r)-CDW phase is obtained 
for V±, J± intermediate, whereas for strong J±, the (tt, 0)-BOW is obtained and for strong V±, the (|, ^)-CDW is 
obtained. The phase diagram for Kp~ — 0.3 is sketched on Fig. 1141 
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FIG. 12: The topology of the phase diagram of the quarter filled ladder for 1/2 > Kp- > 1/3. The (f , f )-CDW appears at 

large V±. 



D. Commensurate-Incommensurate transitions 

1. Deviations from quarter filling 



Away from quarter filling, the number of particles is fixed via a chemical potential that affects only Hp+. The 
charge modes are described by the Hamiltonian: 



/dx 
2^ 



«,+i^,+(7rn,+)2 + ^(a 



2 VI a 
(2^ 



j dx cos A(j)p+ (92) 
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FIG. 13: The topology of the phase diagram of the quarter filled ladder for Kp- < 1/3. The size of the (|-,7r)-CDW phase 
region has been exaggerated. 
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0.01 0.02 0.03 0.04 



FIG. 14: The topology of the phase diagram of the quarter filled ladder for Kp- < 1/3. The (tt, 0)-BOW is obtained for strong 
J± and the ^)-CDW is obtained for strong V±. The (7r,7r)-CDW is obtained in the intermediate regime. 
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where fj, measures the difference in chemical potential A standard argument22ii2Sii2i then shows that for a difference 
in chemical potential of the order of the gap of the quarter filled system, a commensurate-incommensurate (C-I) 
transition occurs and the system becomes gapless. The renormalized K*_^_ exponent at the gapless point is suchifii 
that the dimension of the operator cos40p4_ is one, which gives K*_^ — 1/A. This result has been previously derived 
from different arguments in Ref. and checked against numerical simulations. 

In the incommensurate phase, the |-CDW correlations are the dominant ones as they behave as (e*'^p+(^)e~*'^f+("^) ^ 
x'^p+Z"^. Subdominant vr-BOW or tt-CDW correlations are also present and behave as: (e'^'^''+(^)e~*^'^''+^°'') ^ 
x^'^^p+. Right at the transition, the respective exponents are 1/8 and 1/2. 

In principle, superconducting correlations are also possible. The only surviving superconducting order parameter 
is the d-wave one. Its lattice expression is given by: 



Oici'i') = ^{Ci,l,aCi,2,-a " C^^2,crCz,l-a)- (93) 

(T 

One finds that this pairing operator behaves as e'^f+''^''e~'^f+'*'^ ~ 2;-i/(2^fp+)^ Right at the C-I transition, these 
correlations have an exponent 2 which indicates that superconductivity is largely dominated by CDW correlations. 
In a doped quarter filled insulator the presence of the spin gap is therefore insufficient to render superconducting 
correlations dominant. This is to be contrasted with the case of the half-filled ladderiiS^ where superconducting 
correlations are dominant in the conducting phase. 



2. Effect of a magnetic field 

An applied magnetic field couples to the ladder via a term: 

h f 

i?mag. = - - / dxdx4>cr+ (94) 

It is well known that a magnetic field applied to a spin gap system can produce a Luttinger liquid like phas»i22i when 
the field becomes larger than the gap. The effect of the magnetic field is different in the case of strong repulsion 
{{(f>p-) = f ) and in the case of weaker repulsion {{4>p-) = f ). For strong repulsion the spin-gap is caused by the 
terms gi or 32- The effect of the applied field is to render these two terms irrelevant. As a result, both (j3a-+ and 
(f>a- have gapless excitations leading to a C0S2 phase in the notations of Reflsil This C0S2 phase contains in chain 
2kF charge density and spin-density wave power-law correlations. One finds in the case of C0S2 phase obtained by 



applying a magnetic field to the (tt, tt)— CDW or the (7r,7r)— BOW: 

{OcDW2k^Jx)OcDW2^^, piO)) - cos(^) cos(ma;) ' , (95) 

{S;{x)S;m ^ +cos(|^)cos(m:r)^;^, (96) 

\riT/^\i ^ cos(mx) ,7rx. Const. 

(^P (^)^P (0)> - Const.— ^ + cos( — )-^, (97) 

where exponents have been obtained from the transformation of Rcf.*^ and Sec. IV A 21 by requiring that at the 

transition the dimension of the relevant operators leading to the spin gap in Eq. H71(l be of dimension onc^°i. In the 
case of the C0S2 phase obtained by applying a magnetic field to the (^, ^)-CDWs, one finds insteadiS^: 

{OcDW2k^Jx)OcDW2kj,,A^)) ~ cos( — )cos(mx)— (98) 

{S;{x)S;m ^ 4 + cos(ma;)^^ (99) 

trt-r, \r,Trr.\\ ^ COs(mx) , TTX , CoUSt. 

{Sl{x)S;(0)) ^ Const. \ > + cos( — ) —. (100) 



For weaker repulsion, the spin gap is caused by the terms g4:,g5,g6- The application of the magnetic field then 
produces only a suppression of the gap in giving a COSl phase. The behavior of the induced spin density wave 
correlations depends on whether (jj^- or 6^- is ordered in the parent case. In the first case, which corresponds to the 
(■|,7r) CDWs under strong magnetic field, critical correlations develop in Sp{x) as well as in the {^,qy)-CUW order 
parameters with an exponent of 1. In the second case, which corresponds to the (tt, 0)-BOW and the (tt, 0)-CDW 
under strong magnetic field, only S'f''*' become critical with an exponent of 1/4. 
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VII. CONCLUSIONS 

In the present paper, we have studied charge ordering in the two-leg Hubbard ladder at quarter filling. Focusing 
in the regime of strong-coupling on-site Coulomb repulsion [/, we have investigated the interplay of the interchain 
Coulomb repulsion V± and the exchange interaction J± on charge ordering, and a variety of spin-gapped charge 
density waves and bond-order waves have been obtained. In particular, when the intrachain repulsion is strong 
enough, the ground state of the system exhibits a zig-zag charge order state similar to the phase described in numerical 
studies^. We have obtained the complete phase diagram in the J± — V± plane by numerical integration of perturbative 
renormalization group equations and discussed the transitions between the various charge ordered and bond-ordered 
states. The results show that phase transitions can occur by an ordering in the antisymmetric charge sector or the 
spin sectors. The quantum phase transition in the spin sector, as in the half-filled case, is described by the 0(3) Gross- 
Neveu model^i^ with a mass term and can be either second or first order. The transition in the antisymmetric charge 
sector which is proper to the quarter-filled ladder belongs to the Ising universality class. This type of Ising transition 
is expected to separate the zig-zag charge ordered state from a bond ordered wave phase. We have further analyzed 
the charge and spin excitations in the various gapped phases. Due to charge-spin separation, all of the insulating 
phases have spinless holon excitations of charge ±e. However, the magnetic excitations depend on the nature of the 
phases considered. For strong intrachain repulsion, they can be either massive spinous confined in each chain in the 
case of the (^, ^)-CDWs or domain walls of a dimerized effective zig-zag ladder in the case of a zig-zag charge ordered 
state in agreement with Ref. In the case of weaker intrachain repulsion, we have discussed the analogy of the 
excitations spin/antisymmetric charge sector of the quarter filled ladder with those of a half- filled Hubbard-Kondo- 
Heisenberg chain. In the framework of bosonization we have discussed the connection with other effective models of 
the quarter-filled ladder, such as the spin-orbital model s??'?^ . We have briefly discussed the physics away from quarter 
filling, where commensurate-incommensurate transitions can occur. The analysis of correlations functions show that 
CDW correlations largely dominate superconducting fluctuations at odds with half-filled laddersiSS. Finally, we have 
discussed the effect of magnetic fields strong enough to lift the spin gaps and show that the induced spin density wave 
correlations sharply distinguish the different charge ordered and bond ordered phases. 

We would like to point out that although our results do not directly apply to the NaV205 compound, since in this 
material t± = 2t|| , our model is able to reproduce a zig-zag charge ordered state with spin gap. It would be interesting 
to compare the features of the zig-zag state we predict with the one obtained in NaV205. In this perspective, it 
would be interesting to study in details a model in which interchain hopping is fully taken into account following 
the approach ofSfiiil. It would also be interesting to extend our analysis to zig-zag ladder models, since experimental 
realizations of these systems are now available^. 
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APPENDIX A: DERIVATION OF THE UMKLAPP TERM FROM THE LATTICE HAMILTONIAN OF 

THE QUARTER-FILLED LADDER 

In this section, we give a derivation of the cos4(/)p_|_ term in the bosonized charge Hamiltonian (|23|l of the quarter 
filled ladder following Ref. To derive this term, one needs to separate the low energy processes which keep all the 
particles near the Fermi points dzkp from the high energy processes that involve transfer of particles near the points 
±3kF- By eliminating the latter high energy processes, one is left with an effective action that involve only the Fermi 
points. The details of the procedure to eliminate the high energy states are exposed in Ref. The processes that 
give rise to the cos4:(j)p+ are represented diagrammatically on figure [TBI 

The corresponding terms in the action read for diagram (a): 

= ^ E E E /'^-i E E /'^^ E E X 

r = R.L fci-fc' Pl-'^l •''^ k2,k' Vl,"-! -JO j.^ J., p3,c3 JO 

91 Pi-"! 92 r'2'''2 93 "s'^Ei 

((^ki+qi,-r,pi,ai^k[-qi,-r,p[,a[^k[,r,p[,a[dki,-r,pi,ai){Tl)(c]^^ X 
,-i-,p3 ,CT3 -93 ,-r,p!, ,o-^ '^'=3 '^'P'a ''^3 '^'P'> .(^a ) ('''3 )) d , ( Al) 
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fc, 3kf. -kf 

^~ ^ \ ^ ^ ^ ^7 



(a) (b) (c) 

FIG. 15: The diagrams that give rise to the cos4(^p+ term. 



for diagram (b): 



91 92 ''2'"2 93 '^3'"3 



('^L+<Z3 ,-r,P3 ,0-3 '^fe!, -93 -r,p;, ,cr^ '^'=3 ''"'^3 ^'^3 '^^3,r,P3 ,0-3 ) ('^3 )) d , (A2) 

for diagram (c): 

- ^ E E E / ^-i E E / ^-^ E E f '^^K^A^^'^i- 

r=R,L ki,k> Pi-'i "^0 fca.^o ''^•"s "^0 1-3, fc' p3.<'3 "'O 

91 f'l'"! 92 P2'''2 93 "3 '"3 

("^L)+q3:-r,P3,CT3^fc3-?3,-I-,p!,,0-^'^'=3.-''>P3.0'3'''=3,I-,P3,0-3)('^3))rfj (A3) 

where we use notations similar to those of Ref. ^3- The differences with Ref. are the foUowing. The most important 
is the apparition of a chain index p = 1,2 for the fermions. Another difference (very minor) is that we write the 
annihilation operator for fermions with momentum close to kp, Cu and da for fermions with momentum close to Skp- 
For fermions with momentum close to —hp and —Skp we write the respective annihilation operators cl and cIl. We 
also have by definition —R = L and — L = R. We have the following expressions for X, Y, Z: 

Ki = ^^a.^a'^p^p'^Vx^ad-h-.v) (A4) 
Ylt = (C/a<5,,-.'-2V^|)5p^p, +Fia(l-5p,pO (A5) 
= (^<5,,_,,-V||)(5p,p,+yia(l-5p,pO (A6) 

The operators dR,dL annihilate states of high energy. In the low energy limit, these states do not appear and 
should therefore be integrated oulji^. This is the meaning of the symbol (. . .)d- These states are integrated out by 
using the following Green's function for the d fermions: 

{dk,r,p,cr{T)dk',r',p',a'{0)) = - ^^^ SkM' Sr,r' Spy 6a,a' (A7) 

where t is the transfer integral of the extended Hubbard model. This integration gives rise to a term that describes 
the collision of four fermions with 8kp Umklapp. This term reads: 

^ = ^ / dxAl[^0^{,pl^^^^^i;L^p,,,J^^^^^^^ + H. c. (A8) 



If all the fermions belong to the same chain(pi = p'^ = P2 = p'2) , one recovers the term cos4-\/2(/)p^p derived in Ref. IT^ 
For the case {pi,ai) = (l,t), {p'i,<^i) = (1,1) (^2,^2) = (2,t), {p'2j'^'2) = (2,1) (and all the 4! = 24 cases equivalent 
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by permutation), the Umklapp term IjASp can be bosonized in the form: 

S^g y da;cos(2(0i,| + + 02,T + 02a)) = 5 J cos40p+ (A9) 

This leads to the term we derived phenomenologically. Finding the expression of g is only a tedious calculation. The 
final result is; 



V±{U^ + 10V±U + Wl - 4V|| Vl)a 



9 = — ?^T7ro^ ^ (AlO) 



APPENDIX B: PHENOMENOLOGICAL SPIN DENSITY 

In this appendix, we give a phenomenological derivation of the 2nkp Fourier components of the spin density that 
generalizes the equations H18|I - H2U|I . Let us begin with S'^(x). If we have a system of fcrmions with both spins up and 
down, we can write the following expression^^ for the spin density: 

Pa{x) = Pq - -dx(t>a + cos2m(0Q - kpx) (Bl) 



Expressing these quantities in terms of (f>p,(j)cr, forming the difference and taking into account the term cosvS^o- in 
the Hamiltonian, we obtain the expression: 

S^{x) = ^-j=dx(f)a + Bm sin V2m(f)a sm{mV2(f>p — 2kpx). (B2) 

7rv2 



Since the Hamiltonian of spin excitations contains a term cosv80cr, new terms will be generated by fusion of this 
operator with S''{x). For m odd, the fusions will generate a term sin v^^cr • For m even, the fusions will generate a 
term dxiia- This leads to the following phenomenological expression for S^{x): 

S'^ix) = \=dx4>cr A2m,z siu 2m{V2(j)p - 2kFx) 

TTV 2 — ' 

^ m 

+ X]^2m+l,zSin^/2(/)<,sin(2m + l)(\/2(/)p - 2kFx). (B3) 

m 

The usual expression of the spin density (|20|l is recovered for A2m,z = and A2m+i.z = for m > 1. For the S'^'^(a;) 
similar expressions can be obtained. Starting from the phenomenological Haldane expansion of the fermion creation 
and annihilation operators^^: 



-i/'crix) ~ ev^^ ' \ (B4) 

m— — C30 

we easily obtain: 

S+[x) ~ e*^^" y^ gi[(m-m')(V20p-2fcjj2;) + (m+m' + l)V20„] ^gt^-j 



m,m' 



We note that m + m' + 1 and m~ ml have different parities. We see easily that by fusion with cos VS^o-, the terms 
in with m + m' + 1 odd will be reduced to cos \/2(j)a while the terms with m + m' + 1 even will be reduced to 1. 
The expression of S'^{x) therefore reduces to: 

S+{X) ^ e^^^«"^A2,„+l,.e^(2m+l)(v^0,-2fc,.) 
m 

+6*^"" cosy20,^^2m,.e*2"(^'^^-2fcf"). (B6) 

m 

We can check that the expressions we have obtained lead to rotational invariant expression of the spin correlations 
for Ka = 1, provided that Am,z = Am,x for all m. 
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Eqs. (jB3|) - (|B6|l can also be derived from a more physical argument. The Ogata-Shiba wavefunction tell us that in 
the limit of [//t 3> 1, the spin excitations can be described as a "squeezed" antiferromagnetic spin chains, the spins 
being carried by the electrons^. The spin density should therefore be described by the following expression: 



(B7) 



where the Xi are the positions of the electrons along the chain, labelled in such way that xi < X2 < ... < x^. The 
vector Sn can be decomposed into a staggered and a uniform component as Sn — Jn + (~)"^ri- We will assume 
that both Jn and n„ are slowly varying at the scale of the average interparticle distance so that we can write: 
Jn — J{x — a;„), fin — n{x — Xn), the functions J{x), n{x) varying smoothly between the points x„. We note that 
the integral from — oo to +00 of J is the total magnetization operator and also the generator of the rotations. Thus, 
we expect the functions J and n to obey the following commutation relations: 



[r{x),j\x')] = ieabcJ'{x)6{x-x'){a^b), 
[r{x),n\x')] = ieabcn'{x)6{x-x'){a^b), 



(B8) 
(B9) 



that coincide with the usual commutation relation of the generator of the rotations and the staggered magnetic 
field^"'*. Moreover, since the spin excitations of a spinful Luttinger liquid are expected to be described by a single 
gapless boson, as those of a spin-1/2 chain, it is natural to identify J to the bosonizcd uniform spin density of a spin 
1/2 chain, and n to the staggered spin density. The corresponding expression reads: 



J+{x) = + ijy - e'^^" cos V2(t)a 



ttV2 



; n ^ sm 



in\/20cr 



Let us now introduce a function (j){x) such that 4>{xn) — n?:. We can rewrite the delta function as: 



d{x -Xn) = Y^ S{<j>(x) - nn)— = < 

n n m 

and we have: e**^^^^ — {~)'^ ■ This allows us to write: 



dx ' 



S{x) = J{x)Ye 



d(t) 

dx 



dx 



(BIO) 



(Bll) 



(B12) 



Since ipix) must be titt each time that there is a particle, we have that 4>{x) = + 0x(a;), from which we easily 

obtain: 4>{x) = Trpo — V^4>p- Using the bosonizcd expressions l|B10p of J and n in terms of 4>a the expressions (|B3p - 
(jB6|l are then easily seen to be equivalent to ljB12|l . Applying the expression ljB12|l to our problem we see that the 
terms coming from the 4fci? component of the spin density are less relevant and read: 



(cos40p+ 



t>p-)Ji ■ J2 



(B13) 



The contribution of these terms to the Hamiltonian (|22(l can thus usually be neglected being less relevant. However, 
in the case of Jl sufficiently large, a gap may be formed in the modes , cr± at higher energy scale than in p+. In 
that case, the terms we have derived can lead to a modification of the coefficient go of the term cos40p-|_ in H23|l and 
a change of the ground state expectation value of (j)p+ from J to if Jj^ is large enough. 
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